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The mechanism that drives a metal-insulator transition in an undoped quasi-one-dimensional 
Mott insulator is examined in the framework of the Hubbard model with two different hoppings 
t± 1 and t±2 between nearest-neighbor chains. By applying an iVx-chain renormalization group 
method at the two-loop level, we show how a metallic state emerges when both tj_i and tj_2 exceed 
critical values. In the metallic phase, the quasiparticle weight becomes finite and develops a strong 
momentum dependence. We discuss the temperature dependence of the resistivity and the impact 
of our theory in the understanding of recent experiments on half-filled molecular conductors. 
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The remarkable properties of strongly correlated sys- 
tems near a metal to Mott insulator (MI) transition 
stand out as one of the richest parts of the physics of 
strongly correlated systems This takes on particular 
importance in low dimensional materials at half-filling, 
and especially in one dimension, where spin and charge 
excitations are well known to be invariably decoupled, 
an effect that receives experimental confirmation in the 
one-dimensional (ID) oxide material SrCuC>2 0. How 
this picture modifies when the hopping of electrons in 
more than one spatial direction progressively grows and 
a higher dimensional metallic ground state becomes pos- 
sible is a key issue that remains poorly understood. This 
problem finds concrete applications in organic molecular 
compounds, which constitute very close realizations of 
ID systems. This is the case notably of (TTM-TTP)I 3 
and (DMTSA)BF4 which, as genuine half-filled band ma- 
terials due to the monovalent anions 1^ and BFJ , see 
their Mott insulating state being gradually suppressed 
under pressure [H, 0] ■ Quite recently, resistivity measure- 
ments on (TTM-TTP)l3 at high pressure revealed that 
the temperature scale for the Mott insulating behavior is 
suppressed by an order of magnitude down to Tmi ~ 20 
K at P = 8 GPa of pressure, with a metallic ground state 
expected to occur above 10 GPa It is the objective 
of this letter to propose a theoretical description of this 
transition. 

Bosonization, renormalization group (RG) approaches 
to the ID Hubbard model in weak coupling and its ex- 
act solution from the Bethe ansatz 0,0, Hi, show that 
electron-electron umklapp-scattering processes are a key 
ingredient that promotes the existence of a Mott insu- 
lating ground state at half-filling. The difficulty under- 
lying the mechanism of the MI transition in the quasi- 
1D case resides in the fact that the decoupling of spin 
and charge excitations on the metallic side of the tran- 
sition at higher dimension does not occur, namely when 
Fermi-liquid quasiparticles excitations appear. This issue 
has been addressed theoretically by an RPA treatment 



of interchain hopping |9( , in which its feedback effect on 
the Mott gap is neglected for the self-energy of the one- 
particle Green's function. The RPA results shows the 
existence of a Fermi-liquid metallic state with electron 
and hole Fcrmi-surface pockets when the interchain hop- 
ping exceeds a critical value. The important role of this 
feedback effect has been pointed out recently from the 
use of dynamical mean-field theory, extended to include 
the influence of ID fluctuations (the chain-DMFT) 
- an approach expected to be workable in the strong- 
coupling regime. In the weak-coupling regime, however, 
the electron-electron umklapp scattering and the char- 
acteristics of its coupling to nesting of the whole Fermi 
surface, play a key role [y, HQ, both for the Mott insu- 
lating state and for the possibility of the MI transition 
that emerges in the quasi- ID case. Transverse disper- 
sion which has a tendency towards the realization of the 
deconfined metallic states @, @, E3] is in general over- 
looked in the calculation of one-particle Green's func- 
tion from mean-field like approaches. The two-loop RG 
approach avoids such an approximation and treats the 
Fermi-surface nesting conditions properly. These can be 
strongly altered in the quasi-lD case, especially in the 
presence of frustration in the electron kinetics. 

In this Letter, we apply an A^-chain RG approach 
to quasi- ID half-filled system as developed recently at 
the two-loop level [HI, [12] ■ The momentum dependence 
for nesting of the whole Fermi surface and that for cou- 
plings, including umklapp scattering 13, 14, [lij], are 
taken into account in a systematic way in the calculation 
of one-particle Green's function and four-point vertices 



111 ]. The combined impact of correlations and nesting 



frustration induced by interchain hopping on the trans- 
verse momentum dependence of the quasiparticle weight 
is obtained in weak-coupling regime. 

We consider the quasi-lD half-filled Hubbard model 
on an anisotropic triangular lattice [Fig. HJa)] , with the 
transfer energies t\\ ^> \t± 1 \, \t± 2 \ is the energy along 
chains and txi and t±2 are those between chains). Our 
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FIG. 1: (Color online), (a) Lattice geometry of the present 
model, (b) The corresponding Fermi surface where the case 
for N± — 16 is shown. 



Hamiltonian, with the on site Coulomb repulsion U is 
given by 

H = -h\ E ( c 5a-s-+m,- + 1 L, \) 

- t±i ( c ii,s c j,e+i,s + H - c 

- t±2 ( c jA« c i+i,<+i,« + 11 
+ u J2 n iM n jAi ~ 4as c jAs- (!) 

3,1 3,ts 

The operator Cj < s denotes the annihilation of an elec- 
tron on the jth site in the £th chain with spin s, and 
n j,e,s = Cj t s Cj g s — |. Here ^ = 1, • • • , iVx is the chain 
index. We take the continuum limit along the chain di- 
rection, whereas in the transverse direction we consider 
a finite system having an even number of chains N± with 
the boundary condition Cj t N ± +i,s — c j,i,a- By applying 
a Fourier transform, the kinetic term can be rewritten 
as H = J2k s e (' c ) c l(k) c s(k), where k = (ku,k±) and 
the energy dispersion is given by e(fe) 



-2t\\ cos fen — 

2t^iCosfc^ — 2t^ 2 cos(fc|| + k±) — j«. For \t±{\ <C tu, the 
right- and left-moving electrons in the ID case are well- 
defined due to an open Fermi surface [Fig. QXb)] and the 
shape of the Fermi surface can be parametrized by the 
transverse momentum k± (llT |. To lowest order in the 
interchain hoppings t±_i and ix2> the Fermi surfaces for 
right (+) and left (— ) moving electrons, as a function of 
kj_, are given by kp(k±) — ±(tt/2) ± (t±i/t«) cosk± — 
(t±2/t\\) sinfex and the chemical potential is \i = O(txi)- 
By considering the weak-interacting case and neglecting 
the k± dependence of the velocity, the linearized dis- 
persion used in the RG method takes the simple form: 
£ p (k) — pv[ku — k p F {k±)\ with v = 2tu and p = ±. 

We follow the formulation of the two-loop RG method 
of Ref. [ll| and first introduce the g-ology coupling con- 
stants [see Eq. (2.8) in Ref. namely the backward 
(<7ix)) forward ((fax), and umklapp (<?3x) scatterings with 
opposite spins, and the forward scattering (<7||) and umk- 
lapp (<73||) scatterings for parallel spins. The coupling 
constants are renormalized differently, developing an ex- 
ternal - transverse - momenta dependence in the ver- 



tex corresponding to a patch-index dependence; that is, 
9v ~ > 9u(q±,k± 1 ,kx2)i w here and k±2 are the trans- 
verse momenta for the right-going fermions, and q± is 
the momentum transfer [11] . The magnitude of the ini- 
tial couplings are given by g\\_ = g%\. = #3x = U and 
9 1 1 — .93 1 1 — 0- The RG equations are derived by scaling 
the bandwidth cutoff A(~ 27rfii) as A; = Ae~ l , where 
I is the scaling parameter [ll|. The explicit forms of 
the two-loop RG equations for all coupling constants in 



the case t± 2 — are given in Ref. We solve here 

these two-loop RG equations numerically for a system 
with N± = 16. For even N±, the solution of the RG 
flows indicates the existence of a finite spin gap in the 
low-energy limit, which is expected to vanish in the infi- 
nite N± limit. The characteristic scale In ± above which 
finite size effect would appear can be roughly estimated 
to be l N± w m[A/|txisin(27r/7Vx)|]- For N ± = 16 and 
t±i/t« =0.1 this gives In ± ~ 5, which is sufficient to ob- 
tain the MI transition without finite size effects for the 
present choice of parameters U/t\\ — 2 or 1.5. 

For the bipartite lattice (tx2 = 0), the system remains 



always insulating even for large t±_x U(. This is due to 
perfect nesting condition for the Fermi surface. In this 
case the interchain hopping txi is relevant and becomes 
large under scaling procedure, and some umklapp cou- 
plings become small. However, a macroscopic number of 
these couplings remains relevant due to perfect nesting at 
vector Q — (ir, ir). In the case of large t±2, these relevant 
umklapp couplings are strongly reduced and all umklapp 
terms remain weak. The charge gap then collapses and 
a metallic state with a Fermi surface emerges as a conse- 
quence of nesting deviations that are introduced by large 
interchain frustration. 

One-particle properties are best studied from the 
quasiparticle weight z^ ± . In the RG method at the two- 
loop level [11|, the fcx dependence of the self-energy is 
taken into account in a non perturbative way. This 
contrasts with the chain-DMFT, which treats the one- 
particle self-energy, as a one-chain fcx-independent quan- 
tity of the one-particle Green's function In the 



RG formalism the transverse momentum dependence 11 1 
is taken into account explicitly, and the quasiparticle 
weight of the one-particle Green's function takes the form 
Zk ± = Zk ± z k ± > where z£ and z^ are the contributions 
coming from the normal and umklapp parts of the scat- 
tering, respectively. Here, one can focus on the umklapp 
contribution z% , since z£ remains finite for I < In , in 
the metallic state. The explicit form of zj£ reads [11| 



= exp 



'2N 2 ^ J dlG i<i±-k±K) J HiL,kj_,k' ± 



(2) 



^Su(q ± ,k ± ,k' ± ) — G c(q ± ,k±,k' ± ) + ^c{TT-q ± + k ± + k' ± ,k ± ,k' ± ) 



Gc{q ± ,k ± ,k' ± )Gc{7T-q ± + k ± + k' ± ,k ± ,k' ± ), (3) 



3 



where G c ( q±jk±ik >j = g3X(q ± ,k ± ,ir-k' x )/(2~KV). The quan- 
tity Ju q± k± k * ) is a cutoff function depending on the 
transverse dispersion [ll[ given by J^ q± k± k > ) = 1 for 



! + t 2 ±2 sin 2a. The phase 
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,k' ± ) 

k , = ^o^i[cos/c^+cos(fc^— qx) — cos fc^ 
- 2t [sin A:^ — sm(k± — q±) + sinfc^ - 
The k± dependence of the umklapp- 
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scattering contribution to z k± for several values of txi 
at fixed t_Li/tii = 0.1 is shown in Fig. [2] (a). For weak 
frustration tx2(< t c j_ 2 ), this quantity is small showing 
the existence of an insulating phase with no low-energy 
quasiparticles. On the other hand, for strong frustra- 
tion tx%(> t c ±2 ), it takes sizable values and shows a 
strong fcj_ dependence. In the proximity of the insulat- 
ing phase, the quantity z k presents a broad maximum 
around k±_ « ±7r/2, which behavior implies the emer- 
gence of sections of Fermi surfa ce p , or "cold" regions 
around k w {k p F (±tt/2), ±tt/2) [13]. The dips or "hot 
spots" in the quasiparticle weight correspond to regions 
of the Fermi surface where the nesting, within 2tt/N± 
accuracy, remain the most favorable. In the perfectly 
nested case no Fermi surface is found for any value of 
U. The overall profile is then intrinsically linked to the 
momentum de pen dent nesting properties over the whole 
Fermi surface [15| . These are not taken into account in 
previous analyses 0, [l(J, which reach different conclu- 
sions in the weak-coupling regime. 

The boundary for the metal-insulator transition can 
be determined analytically by noting that the metallic 
phase of the RG analysis is linked to the irrelevance 
of umklapp scattering. The metallic phase boundary 
can thus be obtained when the energy scale of imper- 
fect nesting for umklapp scatterings becomes compara- 



ble to the energy scale 



/tp7exp(-27rf/[/)] of 



the ID Mott gap. By noting that the nesting vector 
of the particle-hole loop, which couples to umklapp, is 
given by Q = (7r, tt ± 2a) with tana = txz/txi, the de- 
gree of the imperfect nesting (the amplitude of the quan- 
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FIG. 2: (Color online), (a) Umklapp-scattering contribution 
of the quasiparticle weight z k± for U/t\\ = 2 and txi/t\\ = 
0.1, with several t±2- (b) Ground-state phase diagram on 
the plane of txi/tu and txz/tu. The dotted line denotes the 
condition of Eq. @. 



tit y K±2a,k ± ,k'J becomes 0= 
boundary for the MI transition is then determined by the 
condition 

t±lt±2 



= cA 



ID 



4- t A 
1 ' ''±2 



(4) 



where c is a numerical constant being of the order of 
unity. The small difference between the numerical re- 
sults and the above analytical expression comes from the 
renormalization of interchain hopping due to the normal 
scattering processes. The Mott gap is also renormalized 
by interchain hopping. The ground-state phase diagram 
in the (£_Li/£i|,tx2A||) plane is shown in Fig. [DJb), where 
the dotted line denotes Eq. (f4]). The ambiguities in the 



cutoff functions of the RG [111 . Il4 ] , prevent us from ob- 
taining a precise location of the phase boundary. 

The Mott transition has also been addressed in two- 
dimensional Hubbard model on an anisotropic triangular 
lattice with nearest- neig hbor hopping t and next-nearest 
neighbor hopping t' 16]. The present model ([1]) can be 
connected to the two-dimensional Hubbard model by tak- 
ing t\\ — > t, txi t, and tx2 — * t' . While our approach 
is restricted to the small interchain hopping, the metal- 
insulator transition obtained here for finite frustration is 
consistent with the numerical results in two dimensions. 

From the solution of the scaling flows of the umklapp 
scatterings, the temperature dependence of the resistiv- 
ity can be qualitatively calculated from the memory func- 
tion approach combined with the RG method by using 
I = ln(A/T) 0, 17 1 - By extending the approach to the 
quasi- ID case, the perturbative expression of the conduc- 
tivity reads p(T) oc N^ 3 E q± , k± ,k> ± & <q± ,k xJ^W e ~ K 
While this formula is not valid for large txi /T [17J] , it will 
depict the qualitative temperature dependence of resis- 
tivity. Typical behaviors of the resistivity for Aj D < txi, 
obtained from this formula, are shown in Fig. [3^a) . At 
high temperature (T > txi, tx2), the effects of inter- 
chain kinetics are masked due to the thermal fluctuations 
where the Tomonaga-Luttinger (TL) liquid behavior is 
reproduced. At low temperature, the insulating behav- 
ior can be seen for small tx2, while the metallic behav- 
ior is found for strong frustration (large txi)- A finite- 
temperature phase diagram, as a function of tx2, can be 
obtained from this behavior, as shown schematically in 
Fig. Elb). At high temperature, the TL liquid state is 
realized, which is followed in the crossover region by the 
development of a Fermi surface. The effect of the frus- 
tration is yet masked by thermal fluctuations, i.e., this 
region can be described effectively by a nested Fermi sur- 
face. At further low temperature, the state moves to the 
Mott insulator or the Fermi liquid (FL) state depending 
on tx2- In the FL state, the effect of tx2 becomes promi- 
nent and the nesting conditions of the full Fermi surface 
are altered. It can also be seen from Fig. OHa) that the 
characteristic temperature Tfl at which the crossover to 
the metallic FL state takes place increases with tx2- For 
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FIG. 3: (Color online). (a)Temperature dependence of the re- 
sistivity for U/tu — 1.5 and t±i/t« = 0.1, with fixed t±2/t» = 
0, 0.02. 0.04, 0.06, 0.08, 0.10 (from top to bottom), and 
po = p(A). T p is the characteristic temperature of the Mott 
insulator, at which the resistivity takes a minimum. Tfl is 
the crossover temperature to the metallic Fermi liquid, where 
the power of p(T) exhibits a change (guided by the eyes), (b) 
Schematic illustration of the t±2-T phase diagram. Tp L is the 
bare temperature scale for imperfect nesting [left of Eq. Q], 
which is reduced to Tfl due to the correlation effects. 



Ap D > t±x, on the other hand, a direct crossover from 
the TL liquid state to the Mott insulator would be seen, 
and neither the FL state nor the crossover region appears 
in the phase diagram. 

We would like to discuss here the impact of our re- 
sults on the understanding of the phase diagram of 
the quasi-lD molecular compound (TTM-TTP) K un der 
pressure The extended Hiickel calculations jlij in- 
dicate that there are two kinds of interchain transfer in- 
tegrals, namely t±i « 9 meV, and t±2 ~ 6 meV, which 
are small compared to t|| « 260 meV. This emphasizes 
the pronounced hopping frustration of this quasi-lD com- 
pound. The Coulomb repulsion between electrons on 
the same molecular orbital of TTM-TTP is weak, and 
the estimated magnitude, U ~ 0.57 eV leads to the 
magnitude U /t\\ k, 2.2 and then the small ID Mott gap 
Ap D /<|| « 0.24. From these figures, the energy scale for 
imperfect nesting [left-hand side of Eq. (|4])] turns out to 
be about 5 meV. As for the magnitude of the charge gap 
at ambient pressure, we obtain A* k, 60 meV, in fair 
agreement with the measured activation energy of resis- 
tivity at ambient pressure The bandwidth goes up 
under pressure, which decreases the ratio U /t\\ and in 
turn the charge gap. To reduce the gap down to the 
scale of imperfect nesting, one has roughly to double 
the hopping amplitudes. This represents a reasonable 
increase of band parameters under 8 GPa of pressure, 
and is consistent with the existence of a MI transition 
in (TTM-TTP)I 3 Q. Finally, we note that at ambient 
pressure, below T c ~ 120 K, this compound develops a 
nonmagnetic state (likely of the spin-Peierls type), which 
comes with an intra-molecular charge disproportionation 
However, since the resistivity already shows an insu- 



lating behavior above T c , the present system can be in- 
deed considered as a Mott insulator rather than a charge- 
ordered insulator, and the charge disproportionation can 
be accounted for by dimerization of the tilted TTM-TTP 
molecules. 

In summary, we have examined the effect of interchain 
frustration on the half-filled quasi- ID Hubbard chains by 
applying an 7Vj_-chain two-loop RG method. The tri- 
angular lattice geometry of the system is found to be a 
key factor in the stability of the Mott insulating state and 
whenever the alteration of nesting conditions due to frus- 
tration in the transverse hoppings reaches some thresh- 
old a metallic state is restored. Our results find direct 
application to the description of the frustrated quasi-lD 
half-filled compounds. 
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